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Abstract. A formalism for semi-adiabatic cranking mass parameters is presented. 
For the fission process of 234 U, the time-dependent pairing equations of motion 
were used to calculate the excitation energy and to extract values of the cranking 
inertia. A fission barrier is determined by minimizing the action trajectory in a 
five dimensional configuration space spanned by elongation, necking, deformations of 
fragments and mass-asymmetry. The deformation energy is computed in the the frame 
of the microscopic-macroscopic model. The two center shell model with Woods-Saxon 
potentials is used in this context. Values of the inertia for excited fissioning systems 
are reported. A dependence between the cranking mass parameters and the intrinsic 
excitation energy is evidenced. 
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1. Introduction 

In microscopic-macroscopic treatments of nuclear fission PQ, the whole nuclear system 
is characterized by some collective coordinates associated with some degrees of freedom 
that determine approximately the behavior of many intrinsic variables. The basic 
ingredients in such an analysis is a shape parametrization that depends on several 
macroscopic degrees of freedom. The generalized coordinates of deformation vary in 
time leading to a split of the nuclear system in two separated fragments. Thus, these 
coordinates describe the change in time of the average field. Traditionally, in the 
large amplitude collective motion like fission or fusion, the calculation of the mass 
parameters are made within the cranking model [21 [3] . These parameters are calculated 
adiabatically. In reality, the nucleus can be excited with respect the ground state at any 
deformation. The excitation energy is shared between the collective kinetic energy of 
the fragments and the intrinsic energy released finally through gamma-ray and neutron 
emission. As mentioned in Ref. pE], no unique separation between these two kind of 
dumping can be realized. Several attempts were realized in order to investigate the 
change of the cranking mass parameters as function of the intrinsic excitation energy. 
In this respect, investigation based on the cranking model have been done for the 
rotational motion [5] and for the behavior of the temperature dependent mass parameter 
P, [71 E] • A theory for non-adiabatic cranking was proposed in Ref. [9] where the time- 
dependent many-body Schrodinger equation is solved and the time dependence of the 
collective motion is determined with the classical Lagrange equation of motion. It is 
worth to mention that in Ref. [10] an expression is derived for the collective kinetic 
energy containing corrections up to the fourth-order in the collective velocities. The 
fourth-corrections to the cranking approximation at hyperdeformations were found very 
large. The authors concluded that the perturbation treatment when applied to fission 
is invalidated. 

In this paper, a new treatment for the cranking model is realized that takes into 
account the intrinsic excitation energy. In the next section, a semi-adiabatic cranking 
approximation is derived for even-even systems to investigate the general features of 
the cranking mass parameters when the nucleus is internally excited. To make the 
problem tractable, it is assumed that the system deforms slowly in time. Therefore, 
the matrix elements of the time-derivative of the wave functions are smalls. Because 
these matrix elements are responsible for quasi-particle excitations, our approximation 
allows to consider the system in a seniority- zero state during fission. In other words, the 
contributions to the inertia originating from seniority-two configurations are considered 
negligible. In section [31 the mass parameters are calculated. To obtain the inertia, 
a least action trajectory for the evolution of the nuclear system is determined. The 
single-particle energies for protons and neutrons were calculated along this path. 
Then, the intrinsic excitation energy of the fissioning system is evaluated within the 
time-dependent pairing equations of motion using different values of the inter-nuclear 
velocities. At each deformation, instantaneous values of the single-particle densities and 
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of the pairing moment components are deduced. The semi-adiabatic effective mass is 
evaluated using these values. Some behavior concerning the dependence of the cranking 
mass parameters versus the intrinsic excitation energy are extracted. In the last section, 
a discussion is made. 

2. Formalism 

The single particle motion of a particle in an average field is governed by the 
Schrodinger equation that includes a Hamiltonian with pairing residual interactions. 
This Hamiltonian depends on some time-dependent collective parameters q(t) = {q v (t)} 
[y = l,...n), such as the internuclear distances between the nascent fragments, the 
mass- asymmetry, the fragment deformations or the necking parameter: 

H(t) = ^k[q(t)](a^a k + a^a k ) - G ^ ajJajfaiO*. (1) 

fc>0 k,i>0 

where e k are single-particle energies, G is an monopole pairing interaction constant and 
a~l denote creation operators. In order to obtain the time dependent pairing equations 
of motion, we shall start from the variational principle taking the following energy 
functional 

d 

C=((p\H- ih— - \N | if) (2) 

and by assuming the many-body state formally expanded as a superposition of time 
dependent BCS seniority-zero and seniority-two adiabatic wave functions 

| <p(t)) = co(t) | (Pbcs) + ^2cji(t)a+af | B cs) (3) 

3,1 

where 

I 0bcs) = Y[{u k + v k ata±) | 0) (4) 

k 

is the seniority-zero Bogoliubov wave function. The excited seniority-two configurations 
are obtained by mean of quasiparticle creation and annihilation operators: 

Uk = u k a k - v k ap, a k = u k a k + v k a k ; (5) 

a£ = u k a k - v* k a k , = u k a^ + v* k a k ; 

In definitions (j2J) and (j3J), Co and Cji are amplitudes of the two kinds of configurations, 
A is the chemical potential, and N is the particle number operator. Because only the 
relative phase between the parameters u k (vacancy amplitudes) and v k (occupation 
amplitudes) matters, in the following u k is considered to be a real quantity and v k a 
complex one. To minimize the functional, the expression §2§ is derived with respect the 
independent variables v k , v k , Cq, Cji (the amplitudes of the wave function), together with 
their complex conjugates, and the resulting equations are set to zero. 

As deduced in Refs. [HI [12] and as detailed in the Appendix, for a seniority- zero 
nuclear system, the pairing equations of motion are: 

ifrpk = K k A* - n* k A ; 

ihk k = (2p k - 1) A + 2n k (e k - A) - 2Gp k n k ] (6) 
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This system is sometimes called time dependent Hartree-Fock-Bogoliubov equations 
[12J. These equations were already used to determine the intrinsic excitation energy 
in fission [T3(, Ej or to investigate the pair-breaking mechanism [15]. The following 
notations are used in Eqs. ([6]) and in the remaining part of the article: 

A o = G n k ] Aji = G ^ n k ] 

k k^j,l 

K k = u k v k ; p k =| v k | 2 ; (7) 

where Ao is the gap parameter for the seniority-zero state while Aji are related to 
seniority-two states. n k are pairing moment components and p k are single particle 
densities. 

An estimate of the intrinsic seniority-zero state excitation energy [11] can be 
obtained with the relation 

E* = E — Ebcs (8) 

where 

E = (0bcs \H-\N\ 0bcs) = 2$> fc (e fc - A) - i-^-L (9) 

k ^ k 

is the expected value of the Hamiltonian (00) for the seniority-zero state and -Ebcs is the 

stationary energy obtained by replacing n k and p k within the time-independent BCS 

parameters k k and p k in formula ([9]). 

As indicated in the Appendix, the next semi-adiabatic cranking formula can be 

obtained for the effective mass parameters B: 



B Ufl — Bi Ufl + B2un 



I Tp Tp \ I K m ^p^\n n \ K n ^p^\K,m\ |2 /™ I OH I „\/„ I 8H I \ 

9fi2 ( E mn-E )\ |>tm|v _ - |Bn|v _ | ( m \ — \n}(n\^\m} 



[E mn — J2k^m,n T k — Eq + J^k T k ) 2 (e m — £ n )' 



( E — E ) ( Km ® pm — Km ® K * m \( K * m — K * m ® Km ) 

+2^, 2 — — ° 9q v ftjjj dq v )^p m 9<Jfj Km <9<} M ' 

m (E mm — J2k^m Tk + T m — E + J2k Tk) 2 

where the values of p k and K k are solutions of the time dependent pairing equations 
(J6j) i and the index v is associated to the collective coordinate q u .. Eji are exactly the 
expected values of the Hamiltonian (JTJ) for seniority-two configurations: 

E jt = («+« r + 0Bcs \ H-\N\ at a +0 BCS ) (11) 

= 2 J2 Pk(e k - A) - - G £ P*+ I e, - A | + | e l - A |; 

^ = (a/ a +0 BCS \H-\N\ a+«t0 BC s) = 2 £ Pk (e k - A) (12) 

I A- I 2 

+ 2u 2 \ej-X \ -—^r- + *Aij + S* A ii " G E I «i I I p)\ 
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and T k are energy terms associated to single-particle states: 

T k = 2 Pk (e k - A) - 2Gp\ + ^ ± "* A ° - l) ; (13) 

The mass parameters can be obtained from formula (11 01) only if the derivatives 
dprn/dq^ and dn m /dq v are supplied. In the BCS stationary state, the derivatives 
dpm/dqu and dk m /dq u depend only on the collective variables q v and the cranking 
formula can be easily obtained. On the other hand, if the system follows a trajectory 
in the configuration space and the motion is characterized by some given collective 
velocities, then the derivatives depend also on q v . Therefore, the second term in Eq. 
f fTUj) depends on the history of the system and can be calculated only if the variations 
of p k and K k are known. The inertia along the trajectory is 



9q v <% 2 ^ {E mm - E 



Krn &Pm Km d^m, 

pm dR K, m 9R 



dR dR m {Emm — J2k^ m T k + T m — E + J2k Tk) 2 

where one coordinate R is taken as the independent variable. The other coordinates 
q v are taken to be function of R [5]. From definition, the collective kinetic energy is 
E c = BR 2 /2 and the last term in Eq. (|T4l) becomes: 

l TP TP \ I K m A 1 2 

o x-^ y^mm ~~ ^0) ~Pm ~ TT^m 
£1 _ ^2 ^ ' ' pm' K m (15) 

m {Emm ~ Hk^m Tk + T m — Eq + J2k Tk) 2 

This term depends only on the derivatives with respect to time k m and p m . Their 
expressions are given by the coupled channel system of equations (J6]). That means, 
the term (j!5p is practically independent on the collective velocity and depends only 
on the values of e m , K m and p m , that is on the structure of the system and its intrinsic 
excitation and represents the minimal collective kinetic energy. Therefore, the term (fT5"l) 
can be viewed as an ground collective kinetic energy. In other words, if the nucleus is 
internally excited and the derivatives of the probabilities k, p are different from zero, the 
system possesses a minimal collective kinetic energy. The term ( JT5l) can be considered 
as responsible for energy transfer between intrinsic and collective degrees of freedom. 



3. Results 

The calculation addresses the fission of 234 U. As already mentioned, the basic ingredient 
in our analysis is the nuclear shape parametrization. The nuclear shape parametrization 
used in the following is given by two ellipsoids of different sizes smoothly joined by a 
third surface obtained [16] by the rotation of a circle around the axis of symmetry. 
Five degrees of freedom characterize this parametrization: the elongation given by the 
inter-nuclear distance R between the centers of the ellipsoids, the two deformations of 
the nascent fragments characterized by their eccentricities £j = \j\ — b 2 j a 2 (i = 1,2), 
the mass asymmetry given by the ratio between the major semi-axis of the fragments 




Figure 1. Nuclear shape parametrization. 



7] = ai/a 2 , and the necking parameter related to the median curvature C = S/R3 (R3 
being the radius of the intermediate circle and 5 the sign associated to the curvature). 
The meaning of all geometric parameters can be understand by inspecting the Fig. [TJ 

As specified in Ref. [5], first of all, a calculation of the fission trajectory in our five- 
dimensional configuration space, beginning with the ground-state of the system up to the 
exit point of the barrier must be performed. This can be done by minimizing the action 
integral. For this purpose, two ingredients are required: the deformation energy V and 
the tensor of the effective mass. The deformation energy was obtained [lj by summing 
the liquid drop energy with the shell and the pairing corrections. The macroscopic 
energy is obtained in the framework of the Yukawa plus exponential model p2] extended 
for binary systems with different charge densities |18j . The Strutinsky microscopic 
corrections were computed on the basis of the Woods-Saxon superasymmetric two-center 
shell model [16J. This model gives the single particle level diagrams by diagonalizing a 
Woods-Saxon potential, corrected within spin-orbit and Coulomb terms, in the analytic 
eigenvalue basis of the two center semi-symmetric harmonic model [T91I20] . The effective 
mass is computed within the cranking adiabatic approximation as given in [5]. After 
minimization, the dependences between the generalized coordinates q v [y = 1, ...5) in the 
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region comprised between the parent ground state configuration and the exit point of the 
external fission barrier supply the least action trajectory. The ground-state corresponds 
to the lowest deformation energy in the first well. The least action trajectory is obtained 
within a numerical method. Details about the numerical procedure of minimization and 
about the model can be found in Refs. [T4"l [T6l [21] and references therein. Plots of the 
minimal deformation energy surface as function of the necking coordinate C and the 
elongation R are displayed in Figs. [2] (a) and (b). In Fig. [2] (b) the minimal action 
trajectory is also plotted. The resulting 234 U fission barrier is plotted on Fig. [2] (c) as 
function of the distance between the centers of the nascent fragments R. Some nuclear 
shapes obtained along the minimal action trajectory are inserted in the plot. 

The realistic level scheme along the least action trajectory were also obtained within 
the superasymmetric Woods-Saxon two-center shell model. Within the energy diagrams, 
the system (jHJ) is solved having as initial conditions the stationary BCS state in the 
ground state configuration. Three values of the internuclear velocity R were taken into 
account: 1 x 10 4 , 1 x 10 5 and 1 x 10 6 fm/fs, corresponding to a time to penetrate the 
barrier comprised between [10~ 18 , 10~ 20 ] s. In Fig. [3] (a), the intrinsic excitation energy 
given by Rel. (JED is plotted as function of the internuclear distance for the three values 
of R. If the internuclear velocity increases, the excitation energy becomes larger. 

The inertia along the trajectory was calculated within three different approaches: 
the adiabatic cranking model of Ref. [5[ [22} [23], the formula (1141) within stationary 
BCS parameters k and p and the same formula within k and p values given by Eq. ([6]). 
The total effective mass is the sum of partial values obtained for neutron and proton 
subsystems. 

In Fig. [3] (c), the inertia along the minimal action path calculated within the 
adiabatic cranking model is displayed with a thick line. In the ground state of the 
first well, the inertia is very small. Large values of of the cranking inertia are obtained 
around the exit point of the outer fission barrier, close to the scission configuration. 
A value approaching the reduced mass is obtained after the scission point. A similar 
behavior, exhibiting very large values of the inertia around the touching configuration, 
was obtained for light systems in Ref. [21], where a version of the Woods-Saxon two 
center shell model based on the molecular orbital approach is used. In the same plot, a 
thin line gives the inertia obtained within formula (fl4"l) by using the stationary values k 
and p. This semi-adiabatic inertia exhibits a similar shell structure as the adiabatic 
cranking model. The semi-adiabatic model gives in general lower values than the 
adiabatic cranking one, excepting the regions of the ground state and of the scission 
point where the values are larger. 

In Fig. [3] (b) a comparison is made between inertia computed within relation (fl4"l) 
for different values of the internuclear velocity. The same values of the internuclear 
velocities R were taken as those used to determine the excitation energies. The inertia 
along the minimal action trajectory has the largest values for the lowest R taken into 
consideration. In general, the values obtained for R = 1 x 10 4 fm/fs represent an upper 
limit for the magnitude of the inertia. It must be mentioned that for this collective 
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Figure 2. (a) Minimal values of the deformation energy in MeV as function of the 
necking coordinate C and the elongation R for 234 U. (a) Contours of the deformation 
energy in step of 1 MeV. The least action trajectory is superimposed, (c) Potential 
barrier. Some shapes obtained during the fission process together with the values of 
the elongation R are inserted. 



velocity, the calculated intrinsic excitation energy is practically negligible. For higher 
values of R, when the excitation energy becomes important, the magnitude of the inertia 
is considerably lower. The values obtained for R — 1 x 10 5 fm/fs and R — 1 x 10 6 fm/fs 
are situated in the vicinity of those obtained within the stationary values. For the three 
collective velocities, the shell structure resemble. 

4. Discussion 

In Ref. [H], a temperature dependent cranking model that includes a parameter for 
the dissipation was used. It was shown that the increase of the temperature smoothen 
and lower the mass parameter. In the high temperature limit, the cranking results are 
quite close to the irrotational flow results. Our calculations also show that the inertia 
decreases when the excitation energy increases. 

Several approximations were derived in the literature to compute the collective 
inertia. In the frame of the response theory [25J the Hamiltonian is expanded around 
a particular value of the macroscopic coordinate and the first derivative is treated with 
the time dependent perturbation approach. Another method is related to the generator 
coordinate method proposed in Ref. |26j within the Gaussian Overlap Approximation 
(GOA). This method allows to obtain a representation of any operator in the collective 
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Figure 3. (a) Intrinsic excitation energies along the minimal action trajectory as 
function of the elongation R for three values of the internuclear distance velocity R: 
full line R = 1 X 10 4 fm/fs, dashed line R = 1 x 10 5 fm/fs and dot-dashed line 
R= 1 x 10 6 fm/fs. (b) The thick lines represent the inertia B divided by the reduced 
mass [i for the three velocities taken into consideration as function of R. The same line 
types are used as in panel (a) for the values of the velocities. The thin line corresponds 
to the inertia computed within BCS parameters k and p. (c) Comparison between the 
inertia calculated within formula (|14p with stationary values h and p (thin line) and 
the classical cranking formula [5] (thick line). Asymptotically the two inertia reaches 
approximately the reduced mass \x. 



space. In general, the GOA mass is about 2/3 times smaller than the adiabatic 
cranking mass [27J but both quantities exhibits a similar shell structure. Adiabatic 
mass parameters for fission were derived also from the Time-Depedent Hartree-Fock- 
Bogololiubov (TDHFB) theory [28]. A comparison between values of the collective mass 
tensor obtained with three diferent models, i.e., cranking, GOA and adiabatic TDHFB, 
showed that the adiabatic TDHFB mass exhibits more pronounced variations than the 
cranking and the GOA masses. As expected, the GOA gives smaller values than the 
cranking model. Our results show that for stationary values of k and p, the inertia is 
lower than the adiabatic cranking values. 

In Ref. [29J the second + collective energy level was calculated for the rare- 
earth and actinide nuclei using the Bohr-Sommerfeld quantization rule. The results 
overestimates the experimental values in average within a factor 2. This discrepancy 
ca be caused either by the shape of the potential in the ground state or either by a 
systematic too low value of the cranking inertia. In our calculations, in the first well 
has larger values than that obtained within the adiabatic cranking model. 
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In conclusion, a semi-adiabatic formalism based on the time-dependent pairing 
equations was described. The values of the inertia obtained within this model exhibit 
a strong dependence on the intrinsic excitation energy. In general, the semi-adiabatic 
inertia have a similar shell structure as the adiabatic cranking model. 
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6. APPENDIX 

In order to simplify the calculations, the variations of the parameter u k and v k as function 
of the seniority number, called blocking effect, is neglected. A treatment involving the 
blocking effect was realized in Ref. [16] in order to generalize the Landau-Zener effect 
and in Ref. [15] to evidence a new dynamical pair breaking effect. Within the wave 
function ([3]), using the notation (00) and the following identities 

d 

(C O 0BCS | 7£ | C O 0BCS) = CqC'o+ | Cq | 2 J2(u k U k + ^k)] (16) 



d 

( J2 c ji a j~ a t <f>BCS | -WL I °kl a k a r"0BCs) = l C *jl C )i+ I C jl | 2 ZJ (^k + V* k V k )] 

j,i¥=j k,l^k k=£j,l 

d 

+ c )i c ki[{u k Uj + v* k Vj){j | — | k)6 a (17) 

k,i^k 01 

- d - 
+ {u i ui + v* l v i )(i I — I l)S jk \; 

d d 
(cjiajaf 0bcs I g£ I c o 0bcs) = 4c (m^ - v d ui)(J I ^ 10; ( 18 ) 

d , d 

(c o 0bcs I ^ I Cjiajajfaca) = -<%Cji(ujVi - v*ui)(j I ^ I 0; ( 19 ) 

d 

(cjtaj affecs I I c o 0bcs) = -c^c^UjVj - v 3 Uj)\ (20) 
d 

(c o 0bcs I I CjjaJ a+faics) = -<\) r jA' l j''j ~ v *jUj)\ (21) 
the energy functional ([2]) becomes 

C =| c | 2 IE 2 I v k | 2 (e k - A) - G | 5> fc v fc I 2 ~GJ2 I I") 

I k k k ) 

+ E I c ji I 2 E 2 I I 2 ( e k - A)+ | ej - A | + | e k - A | 
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G | U k V k \ 2 ~ G I V k 14 

kft,l kft,l 



+ E 



°jj I" \ E 2 I Vk I 2 ( Cfc - A) + 2m 2 | e - Ajj 



k^j 



-iTi |cqCo+ I c | 2 XI ^ 



k^j 

[v* k v k - v* k v k ) 



k+3 



k+i 



+ E 
+E 



4^+ 1 en r e 



c jj c n + I c 



E 



+ 



E [ C *jl C o( u j v l - + CjiC* (v*Ui - V*Uj)] (j 



0_ 

Of 



I) 



+ E r j.i r " ( ' "j'\i + r J<ij) + 'v. r o( + 



+ E E C jl C mn 

- d 

+ {u n ui + v*vi)(l I — I n)5 jm 



d 



(22) 



where the dot represents the time derivative. In order to obtain the equations of motion, 
the previous expression must be derived with respect the independent variables ant their 
complex conjugates and the resulting expressions must be set to zero. The assumption of 
a slowly varying deformation with respect the time will be used. In these circumstances, 
the probability to find the system in a real or virtual seniority-two state is much lower 
than the probability to find the system in a seniority-zero one. These statements can 



be translated in the following conditions: 



'<< 1 and | Co | 2 ~ 1. In this context, as 



prescribed in Ref. [16j and using the derivatives with respect vi and v*, two equations 
are obtained: 

Pi 



ihv* = 2v* m {e m 



A) -G 



u i u i 
2 Ul 



+ (««-^)E«*+ 2 ««r 



(23) 



and another for its complex conjugate. The notations ([7]) are used. From Eq. (T23]) and 
its complex conjugate, the system (jSJ) follows and the expressions involving the product 
between v k and v k apearing in equation ( 1221) are determined: 



T, 



ih 



2MQ-A)-G P 2 + #f|-K Z U 2 



K, 



(24) 



Ao Pi 

2' -' ' " ' 2'1/v; h T 2 

The discussion found Refs. [Tlj reveal that the approximation used to obtain Rel. (1231) 
helps us to describe approximately the effect of the residual interaction on dissipation 



Cranking mass parameters for fission 



12 



and we eliminate only the terms related to the collective kinetic energy. The collective 
kinetic energy will be treated separately by solving equations involving different seniority 
states. 

In order to determine the excitations between configurations the equations (!22j) 
must be derived with respect the amplitudes Co, Cq, c*j, Cjj and c*j. The next three 
equations follows: 

d_ 

ot 



c E Q + ihc -c J2T k -ih^2 c iii. u i v *i ~ v j u j) (j ' \ ^7 I = °i (25) 



k ij# 

d 

CjiEji - ihdji - Cji T k ~ ih[c (-UjVi + VjUi)(j \ — \ I) (26) 

d d - 

+ E Cml(UjU m + VjV* m )(m I — I j) + Y Cjn{u n Ul + VlV*) (fl \ — \ I) = 0] 

Cjj(Ejj — Tfe + Tj) — ihdjj + ihc (ujijj — Vjiij) = 0; (27) 

and three equations for their complex conjugates. The notations (jUJ), ffTTT) . (fT2]) . ([TBI) 
and {^^ were used. The homogeneous solutions are: 

c Q (h)(t) = c oW (0)exp ^- 1 -J\e q - E T *)^) I (28) 

Cii(ft)(*) = q,- { (fc)(0)exp [~ " E ) ! ( 29 ) 

c n(h) (*) = c jm (0) exp (- i ^ - ^ r fc + T^dt j • (30) 

If the system deforms slowly, the contributions of the type Cji(j \ d/dt \ I) in Eqs. 
f[25|) . (126]) and (j2"7|) can be neglected and the solutions for seniority-two states given by 
the Lagrange method of variation of constants are: 

— ifi d 

c ji = p ~ F -l. v T ( v i Ul ~ u i Vl ^i I ~fu I ^ C °W ^ 31 ) 

%J — l^k^j,l 1 k — + Z^fc J fc CI 

c j'i = p v T^T^ F -l. ^ T ~ VjUjjcoih) (32) 

~ l^k^j J-k + 'Lj — ^0 + l^k 1 k 

Finally the probabilities to find the system in a seniority-two state are: 

, |2 ft 2 I ^ ~ ^ I" (j_L|j 2 I . 12 .ooN 

I c i« I = « T^r — ^ 7f F _ i I co(h) | (33) 
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I ' 1 2 

I |2 z, 2 I ^ j' ^7 ^ J I I 1 2 / o /i \ 

I c ii I = n 7^ ^ rp , T p j_sr t\2 I C °W I (34) 

K-^jj ~ l^kjtj 1 k + -Lj ~ + l^k^k) 

The total energy of the system is; 

E= | c | 2 E + J2 I | 2 

= (1 -El 9;i T)^) + E I ^T^ii (35) 

hi 3,1 

= E + E I c j7 | 2 (Efl-Eo) 

3,1 

If we consider that the collective kinetic energy is 

E \ B vAv% = E I c 3i I 2 ( E H - E o), (36) 

v,ii £ j,l 

and substituting the time derivative with derivatives with respect the generalized 
coordinates 

(37) 

and considering that | cq | 2 ~ 1 it follows 

{E mn - Eo) | v m u n - u m v n | 2 (m \ ^ | n)(n | ^ | m) 



£^ = 2ft 2 e 



1 2 ^2 y' (^ mm ( Um %™ ^ m ) ( Um d &z v * m j£ ) ^ 

m (E mm — J^k^m Tk + Tm — i?o + J2k Tk) 2 

where B Ufl are the effective mass parameter. 
Using the identities 

, 8 , < i I # I j > 

< i hr + \ >= — 39 
at €j — €i 

the dependences with respect the derivative of the Hamiltonian are evidenced: 



B ni = 2h 2 E 

m,n^m 



(E mn - E Q ) I v m u n - u m v n | 2 (m | fj | n)(n | | m) 

(-E'mn — Sfc^m,n ^fc — -^0 + Xlfc Tfc) 2 (e m — e n ) 2 



+2 ^2 J2 m 9q " dq » 9q K 9q » (40) 

m [Emm — Hk+m ^k + T m — E + J^fe Tfc) 2 

The values of Uk and t> k are solutions of the time dependent pairing equations. Using 
notations (j7|), the expression ffTOl is eventually obtained. 

It must be noticed that by replacing the energies with the approximate values 



[3D] Eji « ^/(ej - A) 2 + A 2 , + J(ei - A) 2 + A 2 , = Ej + E h by neglecting the differences 
between the sums of T^-terms and by using the stationary values v and u in the previous 
expression, it is straightforward to obtain the usual cranking mass parameter for an 
adiabatic BCS state [5]. In this context, the following identities are also needed: 

% = 2 (1 + ^T ); 
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